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Abstract —Explicit derivation of interferences in hexagonal 
wireless networks has been widely considered intractable and 
requires extensive computations with system level simulations. 
In this paper, we fundamentally tackle this problem and ex¬ 
plicitly evaluate the downlink Interference-to-Signal Ratio (ISR) 
for any mobile location m in a hexagonal wireless network, 
whether composed of omni-directional or tri-sectorized sites. The 
explicit formula of ISR is a very convergent series on m and 
involves the use of Gauss hypergeometric and Hurwitz Riemann 
zeta functions. Besides, we establish simple identities that well 
approximate this convergent series and turn out quite useful 
compared to other approximations in literature. The derived 
expression of ISR is easily extended to any frequency reuse 
pattern. Moreover, it is also exploited in the derivation of an 
explicit form of SINR distribution for any arbitrary distribution 
of mobile user locations, reflecting the spatial traffic density in 
the network. Knowing explicitly about interferences and SINR 
distribution is very useful information in capacity and coverage 
planning of wireless cellular networks and particularly for macro¬ 
cells' layer that forms almost a regular point pattern. 

Index Terms —Wireless cellular networks, Interferences, 
Hexagonal model, SINR distribution. Performance analysis, Ran¬ 
dom user locations. 


I. Introduction 
A. Background and Motivation 

I NTERFERENCE in wireless networks is the prime concern 
of telecommunication actors which continuously attempt 
to minimize it in all the stages of the technology conception 
and deployment (from the standardization to the process of 
network design, planning and exploitation). Interferences are 
often attributed to co-channel interferences which are the 
result of the scarce spectrum reuse in different cells of the 
network. In the network design phase, planning engineers try 
to understand the behavior of the system in terms of radio 
capacity and relate it to interferences. Often, they need to have 
a quick answer on its estimation without having recourse to 
a simulator. This answer should come up with a closed form 
expression of interferences, required to be directly used in 
some tools essential for capacity and coverage planning, such 
as link budget tool. 

In literature, co-channel interferences are often represented 
by the Interference-to-Signal Ratio (ISR) metric [1], named 
also interference factor in 3G studies [2], [3]. It is defined 
as the normalized intercell interferences against the received 
useful signal (received power from the serving cell). ISR 
metric constitutes the chestnut of the performance analysis, in 
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terms of SIR (Signal-to-Interference-Ratio) and SINR (Signal- 
to-Interference-plus-Noise-Ratio). Its evaluation in wireless 
networks has been widely addressed but it is still a fresh topic 
interesting both academic and industrial communities [1]—[11]. 
Interferences depend on the geometric placement of radio sites 
in the network. The two most frequently considered geometry 
models are hexagonal, named also equilateral triangular lattice, 
and random networks. The first assumes that sites are placed 
on a regular infinite hexagonal grid whereas the second 
supposes that sites are randomly distributed according to a 
spatial point process, e.g., Poisson Point Process [6], Although 
random models have recently drawn more attention of the 
research community, hexagonal model is always used as the 
basis for coverage dimensioning (estimation of the number of 
sites required to cover an area with a given quality of service) 
and by operational staff for conducting strategic studies of a 
typical regular network assuming a complete sites’ deploy¬ 
ment. Despite its frequent use in network geometry modeling, 
hexagonal model has been widely considered intractable and 
the evaluation of interferences requires extensive numerical 
computations with simulations [12], [13]. This intractability 
was the key driver in the migration to random network 
models. Besides, system performances are strongly linked 
not only to the network geometry but also to user location 
distributions. With random models, it is hard to analytically 
evaluate the co-channel interference in each user location [1], 
only its distribution and eventually its average are determined. 
Moreover, random network models are often taken stationary: 
distributions of ISR and SINR are independent from the user 
location. 

In this work, we confine ourselves to mathematically tackle 
the tractability issue of the hexagonal model considering 
different user location distributions. Besides, this paper finds 
its motivation from the limit of random network models in 
dealing with non-uniform user location distribution that has 
recently come to prominence [14], 

B. Related works 

Several works have been made to give approximations of 
interferences in regular hexagonal networks based on observed 
simulation results [2]-[4], [11]. To the best of our knowl¬ 
edge, no closed exact formula, mathematically and rigorously 
proved, was found. 

By way of examples, Chan et al. provided in [5] an ap¬ 
proximation of the inter-cell interference distribution assuming 
that the traffic follows a Poisson model. In [6], Haenggi and 
Ganti provided a good review of interferences in regular and 



2 


random networks and in particular they gave lower and upper 
bounds of the cumulated interferences in hexagonal lattices 
assuming that the receiver is placed at the origin [6, p. 19]. In 
the search of the best approximation that well fits the results 
of simulations, Karray provided in [3] a comparison between 
some approximations of ISR in regular hexagonal networks. 
To make the hexagonal network model tractable, authors in 
[2] transformed the hexagonal model to the fluid model: 
interfering sites are continuously and uniformly distributed 
in the plane. It was shown that the fluid model is a weak 
approximation of the hexagonal network. Almost all works 
about interference evaluation in fluid model assume that sites 
in the network are omni-directional. 

Unlike in the hexagonal model, interference analysis is 
tractable when sites are organized according to a Homoge¬ 
neous Poisson Point Process (HPPP), i.e., the number of sites 
in a given area follows a Poisson distribution with constant 
sites’ intensity [6], [8]—[ 10], [15]. An excellent mathematical 
theory of interferences in Poisson network of interferes can 
be found in [6], [8], [9], Despite its tractability, HPPP model is 
very irregular and can not always fit with the geometry of real 
wireless networks because, in reality, sites exhibit repulsive 
behavior while limiting the dispersion of each site position 
within a small area having fixed radius. This is related to 
the fact that radio engineer, while designing sites, places the 
site at its theoretical hexagonal position in a planning tool, 
makes coverage prediction and then looks for real candidate 
sites within a fixed radius from its theoretical position. The 
deviation of the real network structure from the hexagonal 
model is therefore linked to the availability of candidate sites 
(within the search radius) that fulfill certain constraints im¬ 
posed for example by engineering rules, terrain imperfections 
and government charters. 

To obtain performances of regular networks, while profiting 
from the tractability of the HPPP model, Haenggi numerically 
showed in [1] and in [16] with Guo that the SIR and SINR dis¬ 
tribution for hexagonal (or even perturbed hexagonal) model 
with uniform user location can be approximated by shifting 
the SIR and SINR distribution of the HPPP model. In the 
same direction, Deng et al. introduced in [17] the Ginibre 
Point Process that respects a minimum distance between sites 
and thus ensures more the repulsion of sites, feature naturally 
identified in real wireless networks. 

C. Contributions 

In this work, we treat fundamentally the problem of in¬ 
terference calculation and give explicit and exact formulas 
of the downlink ISR in hexagonal wireless networks with 
omni-directional sites. We prove that the ISR metric, as a 
function of the mobile location m = re l6 , admits an absolutely 
convergent Fourier series on 0 and an analytic expansion on 
r. The average of the ISR over the angle 6 is rigorously 
derived without any approximation and involves the use of the 
well-known Gauss hypergeometric functions [18] and Hurwitz 
Riemann Zeta functions [19]. Different identities and accurate 
approximations of ISR are also given and compared with 
other approximations in [2], [3], The explicit ISR expression 


in omni-directional networks is used to derive an accurate 
formula of ISR in tri-sectorized networks. It is also extended 
for any frequency reuse pattern. All the provided expressions 
of the ISR are valid for a propagation parameter b > 1. 

The second contribution of this work is the explicit derivation 
of the network performances in terms of SINR distribution 
in hexagonal networks with an arbitrary user locations’ dis¬ 
tribution reflecting the spatial traffic density [20], The finding 
of SINR distribution entails of course the inversion of the 
SINR, taken as a function of location m. To show the utility of 
the explicit formula of the SINR distribution, we numerically 
investigate the SINR distribution for two scenarios of user 
locations: uniform and heavy-tailed distributions (e.g., Log¬ 
normal user location distribution). 

D. Paper organization 

The paper is organized as follows. In Section II, system 
models, notations and preliminary definitions, including the 
location of ornni and tri-sectorized sites, are given. In Section 
III, we derive the explicit formulas of ISR for omni-directional 
network and we generalize it for tri-sectorized one. This 
includes also the presentation of some new approximations 
and the sketch of their validities and comparison with other 
approximations in literature. Considering a general frequency 
pattern and accounting for shadowing effect in the ISR calcu¬ 
lation are also provided at the end of section III. In Section IV, 
analytical expressions of the SINR distributions are given for 
any user location distribution and numerically presented for 
uniform and non-uniform user location models. Conclusions 
summing up important findings are drawn in Section V. 

II. System models, notations and definitions 
A. Location model of omni-directional sites 

We consider a hexagonal wireless network composed of an 
infinite number of radio sites placed on a regular hexagonal 
grid as shown in Fig. 1 left part. In this latter, each site 
is located at the center of the hexagon and composed of 
only one sector (named also cell). Antenna in each site is 
assumed to have an omni-directional radiation pattern (uniform 
antenna gain in all directions) and covers a geographical area, 
named Voronoi cell. A network with such omni-directional 
sites is called regular omni-directional network. As largely 
considered, we assume also that all cells transmit with the 
same power level P. 

We identify R 2 with the complex plane C and we divide it 
into equally 6 regions f li (0 < l < 5), defined by 

{in CC y < arg(m ) < -—(1) 

Interfering sites are arranged into rings of increasing radii. 
Rings of sites, surrounding the central site located at the origin 
of R 2 , are hexagons. We denote by k, k > 1, the index of 
each ring in the network. The number of sites in the ring k 
and belonging to region is k. We give so the notation Spt,j 
to the site indexed by j in ring k and located in region 12/. 
The site located at the origin of R 2 is denoted formally by 
So,o,o t> u L t0 simplify notations, we omit the lower index in 
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Fig. 1: Hexagonal network model with omni-directional sites 
(Left side) and tri-sectorized sites (Right side). 


.So, 0,0 an d use simply S. We shall use the same label Si y k,j to 
refer to the geographical position of the site. 

Lemma 1: In a regular hexagonal network, the site position 
for fc > 1, is given by 

Si, kd = (2) 

where Dkj = S^/k 2 + j 2 — jk, dkj = '</»»( jL ’ ), <5 is the 
inter-site distance and i = \f^l is the imaginary unit of the 
complex plane C. 

The proof of the Lemma is simple using the geometry of the 
network and it is not provided here. 

A mobile location in the plane is denoted by the random vari¬ 
able m. Likewise, label to is used to denote the geographical 
location of the mobile. In the complex plane, m = re l6 ; where 
r = \m\ is the absolute value of m representing its distance 
to the central site and 6 is the angle coordinate of the mobile 
relative to the horizontal axis. 

For hexagonal omni-directional network, we assume that lo¬ 
cation to is connected to central cell S covering a disk of 
radius R < S so that location variable m < R. The choice of 
R is problematic because the disk circumscribing the hexagon 
over-estimates interferences whereas the disk inscribed in the 
hexagon under-estimates interferences and generates coverage 
holes. For all the theoretical study, we take an arbitrary R, such 

that R < 6, but for the numerical results we use R = 
which is the radius of the disk having the same area of the 
hexagon representing the cell. Actually, the real use of R 
appears only in section IV where we provide the distribution 
of SINR assuming that m is a random variable having dt(m ) 
as its probability measure. 

B. Location model of tri-sectorized sites 

Unlike for omni-directional sites, tri-sectorized sites are 
located at the corner of the hexagon and composed of 3 sectors 
covering each one a cell. Sites are arranged to form always 
a hexagonal grid (see Fig. 1). This model is called regular 
tri-sectorized network. Since each site is divided into three 
sectors, we shall use Si y k y j yC to identify the sector c, 1 < c < 3, 
of the site j in the ring k and located in region fl;. The sector 
c in the central site is simply denoted S c . The azimuth of the 
antenna, in which the radiation is at its maximum, is taken 
relative to the geographical North. Let (f>i y k,j = Cirg(jn—Si y k y j) 
be the angle between site Si y k y j and mobile to, then the angle 


between the azimuth of sector c in the site Si.k j and 
the mobile to is given by 

4>i,k,j,c = ^(2c - 3) + fi yk ,j (3) 

We denote by G(.) the antenna mask 1 of the sector (relative to 
its azimuth) and we define the antenna mask of a site, denoted 
G s , by the sum of all antenna masks of sectors belonging to 
the same site 
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G s (4>i y k,j) = ^2 G(fi,k,j,c) (4) 

C= 1 

C. Propagation model 

In this study, we use simplified Okumura-Hata propagation 
model [21], i.e., the path-loss between site j and mobile 
user to is expressed as 

Li,k,j{m) = a \m- Si yk j\ 2b (5) 

where b > 1 is the amplitude loss exponent, (i.e., 2 b is the 
pathloss exponent), a is a constant dependent on the type of 
the environment (indoor, outdoor, rural, urban,...). We omit 
also the lower index in Lo,o,o and simply use the notation L. 
With this pathloss model and considering the antenna radiation 
pattern G, the received power from sector Si y k y j yC at location 
m is ’ where 4>i,kd,e is given as in (3). 

Since we perform system level analysis, fast fading effect 
is taken from link level results. This latter feeds the former 
with link level curve that relates SINR to the throughput or to 
the packet/bit error rate (see for example [12], [13]). When 
it is constructed, link level curve accounts for fast fading 
considering different types of propagation environment and 
user equipment speeds and categories. Also, the use of fading 
effect in the link performances finds its motivation from the 
link adaptation, performed in real networks, in which fast 
fading is naturally taken in the relation between throughput 
and SINR. In fact, mobile user sends periodically the channel 
feedback report in terms of CQI (Channel Quality Indicator 
reflecting the SINR and including the fading effect) to the 
enodeB (taking LTE technology for example), and when it is 
scheduled, the latter chooses the Transport Block size and the 
Modulation and channel Coding Scheme (MCS) according to 
the received CQI. So, the throughput is given knowing the CQI 
(or SINR) and the instantaneous imperfection of the channel 
due, inter alia, to fast fading. 

Moreover, it is known that, for perfect channel (additive White 
Gaussian Noise Channel), the relation between throughput th 
(normalized with the bandwidth) and SINR £ is provided by 
the famous Shannon’s formula th = log 2 (l + £). To account 
for the imperfection of the channel due to fast fading, one can 
modify Shannon’s formula to have th = K\ log 2 (l + KiC), 
with some constants K\ and A '2 calibrated from practical 
systems as stated in [22]. 

Shadowing effect is not extensively investigated in this 
work but we give in section 1II-D the way of how it can be 

'The antenna mask. G{d) < 1, should not be confused with the antenna 
gain. In logarithmic (dB) scale, antenna gain is the sum of the maximum 
antenna gain and the mask. 
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considered in the ISR calculation. Besides, since the analysis 
is carried out for arbitrary user location distribution, random 
variability of the signal, due to shadowing and fast fading, 
can be aggregated in the user location variable by changing 
the distribution of user locations accordingly. 


D. Definition of ISR 

Definition 1: Let m = re 10 be a location in the plane C and 
let Si t k,j , as defined in Lemma 1, be the location of a given 
site in a regular hexagonal network. We denote by fi,k,j(jn) 
the individual ISR received from site Sz,fc,j i n location m. It 
is calculated relatively to the first sector of central site S by 




L(m) G s (tj)i : k,j ) 
Li, k j(m) G{6- f) 


( 6 ) 


Note that for omni-directional site, there is only one sector 
and hence G s = G = 1. 

With the homogeneity hypothesis that all sectors transmit with 
the same power, fi,k,j(m) is the normalized received power 
from site Si^j, against the useful received power at location 
m served by the first sector of central site S. 


Definition 2: The ISR, denoted by f(m,b), is a function of 
location m and propagation parameter b. It is defined as the 
sum of all individual ISRs 


f(m, b) = Y fi,kA m ) ( 7 ) 

o<;<5,fc>r, 

o<j<fc-r 


III. Mathematical analysis and explicit 

DERIVATION OF THE ISR 

In this section, we use mathematical tools to derive the 
Fourier and Maclaurin series expansion of f(m,b). We first 
analyze the ISR for omni-directional and tri-sectorized net¬ 
work and we point out, at the end of the section, the possibility 
to explicitly account for higher frequency reuse pattern and 
shadowing effect in the ISR calculation. 

A. ISR formulas in hexagonal omni-directional networks 

1) Explicit and exact formulas of ISR: 

Using the definition of the individual ISR in (6) and 
substituting Li k,j by its form in (5), we shall establish the 
following mathematical results. 

Proposition 1: For a given location m = re 1 8 connected to 
the central cell S, the individual ISR, defined in (6), is given 
for omni-directional sites by 

r \ 2b +\ n \ F(& + | n |) e *«(s-®fc,3—if) 
D^J F(6)F(1 + |n|) 

2-Fi (^>,b+\n\,l + \n\, j ( 8 ) 

where 2 -Fi(.,.) and F(.) are respectively the Gauss hyper¬ 
geometric and Euler Gamma functions [18, p. 561]. 

Proof: The proof of proposition 1 is in Appendix A. ■ 


fi,k,j(jn) = Y 


Theorem 1: Let /, defined as in (7), be the ISR function of a 
location m = re 10 and a propagation parameter b in a regular 
hexagonal network with infinite number of omni-directional 
sites. Let 8 be the inter-site distance. Take x = ?, then for 
b > 1 and x < 1, / admits an absolutely convergent Fourier 
series expansion on 6 and an analytic expansion on x 


f(m , b) = Hq(x, b) + 2 'Y / H n (x, b)cos(6n9 ) (9) 


n—1 


where 


<10) 


H n (x,b) 


6F(> + 6 n) x 2b+en 
r(6)F(l + 6n) (1 — x 2 ) b 


and 


U} (b) = 3- b ab)[c(b,-)-ab 1 -) 


(ii) 


( 12 ) 


with £(.) and £(.,.) are respectively the Riemann Zeta and 
Hurwitz Riemann Zeta functions [19, p. 1036], 

To prove theorem 1, we first establish the following lemma, 
proved in Appendix B. This Lemma explicitly provides the 
normalized interferences received at the origin of C in terms 
of known functions. 

Lemma 2: Let z be any complex number such that Il(z) > 
1, where 5ft(,z) is the real part of the complex number z, then 


+oo k— 1 i 

w(z)= §§(* 2 +j 2 -i*r 

=3-*C(s) - C(*,§)) (13) 


where £(.) and £(.,.) are respectively the Riemann Zeta and 
Hurwitz Riemann Zeta functions. 

Proof of theorem 1: 

Using the definition of f(m, b) in (7) and the explicit form of 
frk.jim) in (8), we obtain 


f(m, b)= Y 
0< 

0 

F(6+N) 


0<Z<5,fc>l, n——oo 
0<j<k-l 


r 

D k,i 


26+|n| 


r(6)r(i + |n|) 


2 F 1 \b,b+ \n\, 1 + \n\, 


e in(8-9 ktj -l^) x 


(14) 


Dk, J 


The sum over l is simple to evaluate because it vanishes 
for every integer n not multiple of 6. In addition, the infinite 
hexagonal network is symmetric with respect to the real axis, 
then f(m, b) stays unchanged when we substitute the location 
m by its complex conjugate. It follows that / is an even 
function on 9 and can be written in a closed form 


+OO 

f(m, b) = Ho(x, b) + 2 Y^ H n {x, b)cos(6n9 ) (15) 

n=l 
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where 


H n (x,b) 


6T(b + 6 n) 

r(6)r(l + 6 n) 


+oo k— 1 , \ 2 b+6n 

EE rf: cos(6n9kj) 

k=lj=0 ' k ’ j ' 


2 -Fl 


^ 6,6 + 6 n, 1 + 6 n, 



(16) 


To prove that the Fourier series expansion of / is absolutely 
convergent, we shall explicitly evaluate H n and thus show (10) 
and ( 11 ). 

We first substitute Dkj with 8\/k 2 + j 2 — jk and r with x8. 
Expanding then Gauss hypergeometric function 2 F\ in (16) 
and using lemma 2 yield 


H 0 (x,b) 

and for n > 1 


6x 2b ^ 

r(6) 2 ^ 

v ' h =0 


r(6 + h ) 2 
r(fi + 1 ) 2 


w (6 + h)x 2h 


(17) 


H n (x,b) 


6x 2b+6n y? r(6 + ft)r(6 + 6 n + h) 2h w 

r( 6) 2 r(h + l)T(l + 6n + h) X X 


+oo fe— 1 

EE 


fc=lj=0 


cos(6n8k,j ) 

(fc 2 + j 2 — jk ) b+h+3n 


(18) 


Now when 6 n is sufficiently high, we have the following 
equivalences [18, p. 262]: 


and 


r(6 + 6 n + /i) F(6 + 6 n) . 6 _ 1 

r(i + 6n + h) ~ foTfe) ~ (6n) 


(19) 


+ 00 fc— 1 

EE 


/c=l J=0 


cos(6ndk,j) 

(k 2 + j' 2 — jk) b+h+3n 


~ 1 


( 20 ) 


Replacing the previous relations (19) and (20) in the expres¬ 
sion (18) of H n (x,b) yields 


H n {x,b) ~ 


6 F (6 + 6 n) 2b+6n ST' 2 h 

r(6) 2 r(i + 6n) ^ r(h + 1 ) 


( 21 ) 


The last sum in (21) converges uniformly in the unit disk to 
(t-r'a)!. • obtain then 


H n (x, b) ~ 


6T(b + 6 n) x 2b+6n 
r(6)r(i + 6 n) (1 — x 2 ) b ’ 


as 6 n — > +oo 

( 22 ) 


It follows that the Fourier series expansion (9) of the ISR / 
converges like the series 


H 0 (x,b) + 


I2x 


2b 


E 


r(6 + 6n) 6n 


r(6)(l — x 2 ) b T (1 + 6 n) 


x bn cos(6?i0) (23) 


which converges absolutely for \x\ < 1. In addition, the ISR / 
can be very well approximated by the series in (23) for every 
x < 1 and b > 1 . ■ 


It is interesting to note from theorem 1 the following 
remarks: 

• The explicit formulas of the ISR f(m, b) in (9) and 
its average over 9 in ( 10 ) are rapidly convergent series 


on to , for | to | < 8 , and were unknown before this 
work. The first elements of the series are sufficient to 
capture with a high precision the ISR value at location 
to. Therefore, using these exact expressions is better, in 
terms of computation time, than performing simulations 
with high number of sites. 

• The received interference at location to is l(m,b) = 

f(m, b) —+ 7 . At the origin of R 2 , i.e. to = 0, it is simply 
reduced to the quantity 1 ( 0 , b) = • 

• The Mean Interference-to-Signal Ratio (MISR) for the 
hexagonal model can be simply deduced by averaging 
/(to, b ) over to. When to follows a uniform random 
variable on the disk of radius R = kS (where k could be 

for example 1/2, or l/\/3), the MISR is 


m 


6 V" r(& + h) 2 uj(b + h) 2b+2h (2 a\ 

r ( & ) 2 fro (b + h +l)T{h+l) 2 > 


which is independent from the intersite distance 8. 


Corollary 1: Let /, as defined in (7), be the ISR function 
of location to and parameter b. Let x be as in theorem 1, then 
we have for b > 1 and x < 1 


/(to, b) ~Hq{x, b) - 


12x 


2b 


(1 — X 2 ) b 


2x 


2b 


(1 — X 2 ) b 


E*[( 

;=o L 


1 — xe 


i(6+lTr/3)' 


(25) 


The proof of corollary 1 is simple. We have just to replace 
H n (x,b) in the expression of /(to, b) by the approximation 
of (11), valid for n > 1, and then evaluate the sum in (23). 


Remark: 

• As can be observed from equation (11), H n (x , b) decays 
very fast with n and becomes very small for n > 1. Thus, 
the coefficients H n (x, b), for n > 1, contribute very little 
on the value of /(to, b). The ISR / is then a very slowly 
varying function on 9. This observation was pointed out 
in [3] based on the curve of /. 

• The ISR /(to, b) can be well approximated only with 
H 0 (x,b) with an error equal to 0(x 2b+6 ); where O is 
the Knuth’s notation. 

Corollary 2: For b > 1 and x < 1, the following 
approximations hold. 



f(m , b) « H 0 {x,b) 

(26) 

H 0 (x,b) = 

6x 2b (u(b) + ui{b + 1 )b 2 x 2 + 0(x 4 )) 

(27) 

H 0 (x,b) 

« 6 a ; 26 ( 2 ^ 1 ( 6 , b, l,x 2 ) + ui(b) - l) 

(28) 

H 0 (x,b)z 


(29) 


Proof: Approximation (26) arises directly from the previous 
remark. Approximation (27) results from the Maclaurin series 
expansions of Ho provided in ( 10 ) and we take only the 
development to order 2 with respect to x. 
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Fig. 2: Explicit formula of / in (25) versus simulated ISR at 
each location m. 


To prove (28), let to be the function in (13), take the expansion 
of Ho in (10) and use the approximation u>(b + h) tv 1, for 
h > 1. 

To prove (29), we start from (28) and use Euler transformation 
formula of the Gauss hypergeometric function [18, p. 564]: 

2 F 1 (b,b,l,x 2 ) = (1 - x 2 ) 1 ~ 2b 2 F 1 {l - 6,1 - 6,1, x 2 ) (30) 

To complete the proof, take only the Maclaurin series devel¬ 
opment of order 2 in 2 F\ (1 — 6, 1 — 6, 1 , a; 2 ) because, for x < 1 
and b > 1 , the other terms of higher order are negligible. ■ 


2) Validation with numerical examples: 

In Fig. 2, we plot the explicit formula (25) of the ISR and 
its simulated value in each location m considering different 
amplitude exponent coefficients b G {1.25,1.4,2}; and for two 
cases of angle 6 G {0,7r/6}. For each value of location m, 
the simulated ISR is performed with 1000 rings of cells. We 
notice that the explicit simple form (25) of / is exactly similar 
to the one obtained by simulation following its definition in 
(7), for all x < 1 and b > 1, but in practice we focus only 
on the part where the location m is served by the central cell, 
i.e., x = § < 

It is worth noting that the ISR undergoes small variations with 
9 except at cell edge (x « ^) where its impact becomes 
significant. Moreover, the influence of 6 on the ISR increases 
with b. 

Fig. 3 represents the exact expression of Hq in (10) compared 
with their development of order 0, namely 6u:(b)x 2b , and the 
2 nd order given in (27) for b = 1.1, b = 1.3 and 6 = 2. 
For 6 close to 1, both the order 0 and the 2 nd order match 
well with the exact form of // 0 . However, for high value of 6, 
the exact value of Hq begins to distance from its first orders’ 
developments and mainly from order 0. This is explained by 
the fact that when 6 increases, the coefficient of the 

series expansion of Hq also increases. Hence, we need to go 
for higher order and notably at the cell edge area. To sum up, 
The 2 nd order development (27) is a good approximation of 
H q in the covered area of the cell for 1 < b < 1.5. 


3) Comparison with other approximations in literature: 

In addition to the explicit calculation (10) of Hq, this paper 
gives a simple approximation of Hq in (29) that would match 
very well with Hq for every 6 > 1. To show the effectiveness 
of this approximation, we numerically compare it with other 



-Explicit formula of H Q 

n (10) 

+ Approximation in (29) 


o Approximation with f 

n (31) 

A Approximation with f g 

n (32) 


Fig. 3: Hq and its developments of order 0 and 2 nd order. 


Fig. 4: Comparison of the ISR formula (29) with other 
approximations for b G {1.2,1.4, 2). 


approximations provided in [2] and [3]. The approximation of 
the ISR (provided for the fluid model) given in [2, Eq.15] is 


/,(m> = 73(b^Tj (1 ~ 

Whereas, the provided one in [3, Eq.7] is 


2-2 b 


/2O) = 


((2b — l)x 
(1 — x) 2b 


2b 


1 + 4(1 -x) b + 


(l-x) 2b 
(1 + x) 2b 


(31) 


(32) 


In Fig. 4, numerical results show that, for any value of x < 
the average over 9 of the ISR, Hq, is well captured and 
precisely estimated by the proposed simple approximation 
(29). The approximation f 2 , in [3], is valid for low values of 
x but is limited at the cell edge (with error higher than 20%) 
for all 6 > 1. The approximation /i, in [2], is good only when 
6 approaches 1. For higher 6 (e.g., close to 2), /i moves away 
from Hq and becomes inefficient to estimate the interferences 
in hexagonal network. To reduce the gap between /i and the 
exact value of Hq, Kelif et al. introduced in [2, Eq.16, p. 7] 
a corrective term obtained by simulation. 


B. ISR expression in hexagonal tri-sectorized networks 

Sectorization effect on interferences was largely investigated 
but mostly with simulations [23], [24]. In regular hexago¬ 
nal networks with tri-sectorized sites, interferences obviously 
depend on the antenna radiation patterns of the sectors in 
both serving and interfering sites. In general, both horizontal 
and vertical radiation patterns influence interferences, but here 
we give an approximation of the ISR considering only the 
horizontal pattern. The vertical one is neglected and we assume 
that its effect is taken in the constant term of the antenna gain. 
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As mentioned earlier, antenna pattern (in dB scale) is the sum 
of a constant gain and a mask depending on the angle between 
the mobile location and the antenna’s boresight. 

The antenna mask G s of a site, as defined in (4), is a periodic 
function with period 2 / if all sectors forming the site have the 
same antenna pattern. In tri-sectorized network, the relative 
interference F, received in location to, is the sum over all 
interferences received from each site Si t k,j (for 0 < l < 5, 
k > 1 and 0 < j < k — 1) and from all collocated sectors of 
the same site. 

F{m, b) = -1 + + fs{m , b) (33) 

where the quantity — 1 + corresponds to the intra¬ 

site ISR generated by sectors belonging to the same site and 
f a (m, b ) is the inter-site ISR as defined in sub-section 1I-D. A 
good approximation of F(m,b ) is provided in the following 
proposition. 


Proposition 2: The ISR F(rri. b ), dependent on b and 
received at location to connected to the first sector of site 
S' in a regular hexagonal network with infinite number of tri- 
sectorized sites, has the following approximation: 


F(m, b ) 


G s (6) g 0 f(m,b) 
G{9- f) G( 9- f) 


2a lX 2b ® K^i-se*) 3 ] 
G(B- Vh\l-xe^-m 2b+3 


(34) 


where / is recalled the ISR function for omni-directional 
network, given in theorem 1, and x = ?. The coefficients ao 
and are given by 


a p = — / G s (Q)cos(3pQ)dd , {0,1} (35) 

k Jo 

Proof: The proof is provided in Appendix C. ■ 


It is important to note that the approach, used to give the closed 
formula (34) of the ISR in tri-sectorized network, is applicable 
to any sectorization level such as quadri or hexa-sectorization. 

In order to assess the validity of ISR approximation in tri- 
sectorized network, we use an antenna mask from a real pattern 
having a beamwidth equal to 65°. The coefficients of the site 
mask calculated in (35) for the considered antenna are ao = 
0.62 and ai = —0.19. 

In Fig. 5a, we draw both simulated results and the explicit 
form (34) of the ISR in a hexagonal tri-sectorized network for 
different values of b and for the extreme values of 6: 0 and 

Simulations is performed using 1000 rings of tri-sectorized 
sites. For 9 = 0, r < | whereas for 8 = r varies up to 4/. 
The simple explicit formula (34) approximates very well the 
ISR for all locations in the cell except with small difference 
when to approaches and for high value of b. In addition, 

for 9 = 0, the signal coming from the serving sector is almost 
equal to the signal from the collocated one and consequently 
ISR is higher than 1 even for locations close to site center. 
Furthermore, we observe a high gap between ISR in 9 = 0 
and in # = f • Hence, the impact of 9 is no longer negligible 




0.7 (b) Illustration of the intersite ISR 

G(0-7t/3)f (m) for b=1.4. 


0.3 0.4 0.5 0.6 

x=r/8 

(a) ISR in tri-sectorized network for 0=0 and 0=7t/3. 


o 


Simulated ISR 
Explicit formula in (34) 


b=1.25 


Fig. 5: ISR in tri-sectorized network. 


and it plays a crucial role in the calculation of the ISR. 

In Fig. 5b, the inter-site ISR G{6 — J )f s (m,b ), in (59), 
is represented as a function of 9 and r. In contrast to 
the small impact of 9 in the omni-directional network, we 
notice further that the inter-site ISR in tri-sectorized network 
depends on 9 as much as on r. In fact, the antenna mask 
defines the envelope curve of the ISR calculated in omni¬ 
directional network. Furthermore, Fig. 5b shows that the inter¬ 
site interference is periodic with period This result is easy 
to be mathematically proved. 


C. ISR expression for higher frequency reuse patterns 

Frequency reuse is the core concept of wireless cellular 
networks and consists in using the same frequency in different 
geographical areas, spaced far enough apart from each other. 
To avoid excessive interferences, fixed frequency allocation to 
each cell is performed in the planning process according to a 
frequency reuse pattern, denoted here by v. According to this 
pattern, v different frequency bands of the operator are used 
for each cluster of v adjacent cells. 

Co-channel interfering cells are arranged to form again a 
hexagonal lattice with first ring distant from the cell at the 
origin by A = \Jv8 (see [25, p. 8]), where 8 is the inter-site 
distance and defined as in section II-A. Let /(to, b , v) be the 
ISR of a hexagonal network having a reuse pattern v, it is 
obtained by changing 8 in the ISR of reuse 1 by A. Since the 
ISR of reuse 1, f(m,b) = f(m,b, 1), depends on the ratio 
m/8, it follows that 

/(to, 6, v) = f(m/y/v, b) (36) 

It is clear that f(m,b,v) converges faster when v increases. 
Consequently, for high value of v, the first term of the series 
of f(m/y/v, b) is enough to capture the essential value of 
the ISR. It means that for high value of v, f(m,b,v) ~ 
6ui(b)x 2b /v b , where again x = |m|/<5. 

It is known that for high frequency reuse patterns, the 
operator is obliged to allocate a low bandwidth to each cell 
resulting in a decrease of cell capacity. In order to maximize 
the spatial spectrum reuse while minimizing interferences. 
Fractional Frequency Reuse (FFR) pattern has been proposed 
in GSM networks for the concentric-cells’ deployment (called 
also Underlay-Overlay cells) and in LTE for the Inter-Cell In¬ 
terference Coordination (ICIC) scheme. The fundamental idea 

















of FFR is to apply two frequency reuse patterns in the same 
cell: tight frequency reuse with short reuse distance, restricted 
to users close to the cell, and regular frequency reuse with long 
reuse distance, useful for users at cell edge. For FFR, the ISR 
provided in (36) is still valid but the frequency reuse pattern 
v depends indeed on the location m. A well known example 
of FFR pattern is v(m) = 1(|to| < ro) + 3 x l(|m| > ro), 
where l(|m| < ro) is the indicator function that takes 1 if 
the condition "\m\ < r o" is verified and 0 else and ro is a 
threshold separating cell edge users from those close to the 
cell. This means that a frequency reuse pattern v = 3 is used 
at cell edge whereas at cell center v = 1. 


D. Inclusion of Log-normal shadowing in the ISR calculation 

To take into account Log-normal shadowing effect on the 
calculation of the ISR, the pathloss Lipj(rn) between location 
m and site Sij ; j becomes a Log-normal random variable 
multiplied by its expression in (5). In each location m, 
Li,k,j (m) and (m) are also supposed to be independent 

and identically distributed whenever (. l,k,j) f 
Consequently, the individual ISR, defined in (6), is also a Log¬ 
normal random variable since it is the ratio of two Log-normal 
variables. Let f( sh \m, b) be the ISR when accounting for the 
effect of shadowing, its definition in (7) is changed for omni¬ 
directional network by 


f {sh \ m , b )= E 


126 


0<l<5,k>l, \ m Sl,k,j 
0<j<k-l 


,2b 


Xi,k,j (37) 


where {xi,k,j} is a sequence of independent and identically 
distributed Log-normal variables with mean E[x] and variance 

Var[x\. 

The distribution of the Log-normal sum is not exactly known 
and always constitutes an open problem in probability theory. 
Nevertheless, it has been recognized that the Log-normal sum 
can be well approximated by a new Log-normal random vari¬ 
able using Fenton-Wilkinson method [26]. This is achieved by 
matching the mean and the variance of the new resulted Log¬ 
normal variable with those of the Log-normal sum. The choice 
of Fenton-Wilkinson approach is motivated by its convenience 
for the present study in contrast with other approaches in [27], 
[28], Applying Fenton-Wilkinson method to the Log-normal 
sum in (37), it follows that b) is approximated by a 

Log-normal random variable with mean 

E[fW(m,b)} =f(m,b)E[ X } 


IV. Application to performance analysis: SINR 

DISTRIBUTION 

The performance analysis of wireless networks definitely 
entails the evaluation of the SINR distribution. Its knowledge 
allows to estimate the throughput which is related to the SINR 
by a continuous, differentiable, and strictly monotonically 
increasing function. As mentioned in section 1I-C, this function 
is called the link level curve of the system. 

In this section, we give the explicit formula of the Comple¬ 
mentary Cumulative Distribution Function (CCDF) of SINR 
for a general location distribution. Considering an arbitrary 
user location distribution is of great interest for capacity 
dimensioning because it was shown that the spatial user 
distribution is practically very heterogeneous between cells of 
the network and even inside cells [20], [29]. 


A. Explicit form of SINR distribution 

Before the explicit calculation of the SINR CCDF, we give 
the following definitions. 


Definition 3: The SINR, experienced at location m connected 
to central cell S in a regular omni-directional hexagonal 
network, is the ratio between the useful received power and 
the total received interferences including thermal noise power 

SINR=( V f(m,b) + y 0 x 2b )- 1 (38) 


where yo = ^jf-S 2b , P is recalled the transmitted power of 
the cell, including the antenna gain, a and b are the path loss 
parameters given in (5), <5 is recalled the intersite distance, 
x = |, r) is the average load of the interfering cells and Pjv 
is the thermal noise power calculated in the same spectrum 
bandwidth as for the cell transmit power P. 

Equation (38) is obtained by firstly considering the classical 
definition of SINR at location m 


SINR = 


P/L(m) 

E Vi,k,jP/Lixj(m) + P N 

o,j 


(39) 


where yi,k,j is the percentage of occupied resources (schedul¬ 
ing slots or physical resource blocks) in the interfering cell 
and reflects its load averaged over a period corresponding to 
the averaging window of the measurement reports of mobiles; 
and secondly assuming that all interfering cells have the same 
average load y. 


Definition 4: The CCDF of SINR for any scenario of user 
locations’ distribution is given by the following integral 


and variance 

Var[f {sh) (m,b)] = f(m,2b)Var[x] 

where f(m, b) is the ISR without the inclusion of shadowing 
and given explicitly in (9) or in (25). 

It is important to note that, with the inclusion of shadowing, 
the Voronoi cell associated to the central cell is no longer the 
hexagon or the disk having the same area as the hexagon. 
It becomes the set {m G C | L(m) < Lykj{m),k > 1}, 
instead. 


V{y) = P (SINR > y) = f 1(SINR > y)dt(m) (40) 

Js 

where S is the area of the central cell assumed to be a disk 
of radius R < 5, 1(.) is the indicator function. t(m) is the 
probability measure of the location variable m and reflects the 
spatial traffic distribution in cell S. 

We assume that t is a probability measure on S and conse¬ 
quently 

/ <If(m) = 1 

Js 
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Remark: 

• The CCDF of SINR is widely called coverage probability 
because it gives the percentage of locations in which the 
condition "SINR > y" is guaranteed. 

• When the thermal noise power is neglected with respect 
to interferences, yo ss 0 and the SINR might be similar 
to SIR. 

• When the SINR in each location m is a random variable 
(e.g., through the shadowing effect), definition (40) of the 
SINR CCDF is replaced by 

^f(y) = f P(SINR > y\m)dt(m) 

Js 

In the calculation of the SINR distribution, we consider the 
expression of the ISR in (26) since we showed earlier that, for 
an omni-directional network, the angle 6 has a little impact. 
This latter would be further smoothed when looking for the 
distribution of any transformation of the ISR function. 

Let g be the function defined on [0,1) by 

g(x) = r)H 0 (x) + y 0 x 2b (41) 


where again x = ?, it is clear that g(x) = 1 /SINR. and the 
calculation of SINR CCDF requires the inverse function of g 
which realizes a continuous bijection from [0,p] to [0, </(/?)] 
for every real number p satisfying 0 < p < 1. The inversion 
of g can be explicitly given by series reversion method. In 
the following proposition, we provide an approximation of 
5 _1 (y)- 

Proposition 3: Let g be defined as in (41). Let y = g(x) for 
every x < p = then the following approximation holds 
for x = g~ 1 (y) 


g~\y ) 


C(y,b) 


where 

C(y,b) = ( 


\j\ + \J\ + P(P)C(y,b) 2 

6 bgco(b + 1) 


(42) 


6 yu(b) + y Q 


l 2b and /3(b) = 


6 r)u(b) + yo 


Proof: The proof is given in Appendix D. 


Considering now an arbitrary traffic distribution, we shall state 
the explicit form of SINR CCDF involving the inverse function 
g -1 in the following theorem. 

Theorem 2: Let g be defined as in (41). Recall 5 be the 
inter-site distance and R be the radius of the central cell S. 
Then the SINR CCDF is explicitly given by 

*{y)~T(A(y)), Vy > 0 (43) 


with 


A (y) = min 


Sxg-^-hR 


(44) 


The function T(.) is the marginal Cumulative Distribution 
Function of the location random variable m and is obtained 
from the probability measure t(m ) = t(r,9 ) by dT(r) = 
f Q w dt(r, 9). The function min(u, v) gives the smallest one 
among the real numbers u and v. 


Proof: Taking the definition of SINR and T' in respectively 
(38) and (40) and with the assumption that the impact of 9 on 
the ISR is neglected, we have 

/*2tv nR 1 

^{y)= / > y)dt(r,9) 

Jo Jo 9\$) 

min((5xg _1 (i),J?) r2v 

/ dt(r, 9) 

Jo 
Hv) 

dT(r) 

B. Numerical results for different user locations’ distributions 

In order to evaluate and validate the analytical expression 
of SINR distribution T' for different scenarios of user loca¬ 
tions’ distributions, we consider a uniform and non-uniform 
locations’ distributions in the cell. Note that in practice, the 
estimation of the spatial traffic density (the measure t(m)) 
in each location can be based on processing probes or on 
the manipulation of some key performance indicators of the 
network as in [20]. Furthermore, authors in [29] showed 
with real data measurements that the spatial traffic density 
can be approximated by a Log-normal distribution. Therefore, 
numerical results of this work are provided for both uniform 
and Log-normal 2 user locations’ distribution scenarios. 

Besides, for each location scenario, we simulate 
interferences in a hexagonal omni-directional network 
using Monte-Carlo method and compare it to the theoretical 
results. The considered network is composed of 1000 
rings of cells operating in the band 2600 Mhz of the LTE 
technology. Only 20 Mhz of spectrum bandwidth is available 
and reused everywhere in the network, i.e., frequency 
reuse pattern equals 1. Each cell of the network transmits 
with a power equal to P = 60dBm 3 and distant from its 
adjacent cells by <5 = IKm. We set the radius of the cell to 

R = d\J= 0.5255 and we assume that interfering cells 
are fully loaded, i.e., // = 1. The downlink thermal noise 4 is 
set to —93 dBm. It is worth noting that parameters in dB 
scale are transformed to linear when used for the calculation 
of the simulated metrics ISR and SINR. 

For each scenario of simulation, 20000 users are generated 
according to the considered location distribution in subsections 
IV-B1 and IV-B2. For each generated user location, 
interferences, ISR and SINR are calculated. 

The comparison between theoretical and simulation results 
are provided for 2 types of environment, different by the value 
of the propagation parameter a: Outdoor environment with 
a = 130 dB and deep indoor environment with a = 166 dB. 
Furthermore, We investigate results for 3 different values of 
b: 1.25, 1.5 and 2. 



2 The Log-normal user location distribution should not be confused with 
the Log-normal shadowing. 

3 P = 60 dBm corresponds to 43 dBm from the transmitter power 
amplifier and 17 dB for the antenna gain of the transmitter. 

4 P n = E n + 10 log{W) + NF\ where E n is the energy of the thermal 
noise equal to —174 dBm/Hz, W is the spectrum bandwidth equal to 
20 Mhz and NF is the noise figure of the transmitter, set here to 8 dB. 



10 




°0.5 

DC 

l0.4 


-5 0 5 10 15 20 

SINRin dB 

(a) In outdoor environment 


-5 0 5 10 15 20 

SINRin dB 

(b) In deep indoor environment 


_SINR CCDF with 

Monte Carlo simulations 
o Calculated SINR CCDF in (45) 


Fig. 6: CCDF of the SINR for uniform location distribution 
scenario with different values of b. 


1) SINR distribution for uniform user location variable: 
When the user location distribution is uniform in the disk 
(representing the cell) of radius R, we have dt(m) = r ^i ■ 
The distribution T of SINR in (43) becomes 




A (y) 2 

R 2 


(45) 


Since the ISR is a function of the polar coordinates of m, 
approximating the central Voronoi cell with a disk of radius 
R gives a simple and closed form of SINR distribution. This 
assumption has been widely considered, see for example [30], 
[31]. 

In Fig. 6, we present the curve of 4' in (45) and the SINR 
CCDF obtained by Monte Carlo simulations for both outdoor 
environment (Fig.6a) and deep indoor environment (Fig. 6b). 
The results show that the expression of calculated in (45) 
fits with simulations results of the SINR CCDF quite well. 
Moreover, as expected the impact of 6 on is smoothed 
and thus neglected. In addition, it is important to note, from 
the comparison between Fig. 6a and Fig. 6b, that in the deep 
indoor environment, the impact of the value yo on the SINR 
is clear and tends to decrease it by 2 to 3 dB depending on 
b. This is explained by the high value of the deep indoor 
penetration margin that makes the path loss very high and thus 
the signal arrives very low at indoor locations. Conversely, for 
outdoor environment, yo is low and SINR CCDF should be 
quite similar to SIR CCDF. 

2) SINR distribution for Log-normal user location variable: 
For a Log-normal location distribution, we assume that vari¬ 
able r follows a Log-normal distribution ~ InN(ft, a) with 
mean ft and standard deviation a. Variable 9 is always uniform 
~ U(0,2n). The probability measure of location m is written 
by 

_(i»W-c) 2 , 

, . . e 2 v' 2 dr dO 

dt(m.) = -==——--- 

oV^nQ [“(In(-R) - ft)] r 2 tt 


and consequently the expression (43) of T 1 becomes 


*( 2 /) 


Q [^( ln ( A 0/)) - r)] 

Q[A(ln(J*)-/i)] 


(46) 


where Q is the cumulative distribution function of the standard 
normal distribution and is given by the following integral. 


1 u 2 

Q(z) = _ / e 2 du, \/z G R 

V 27T J — oo 



(a) Traffic hotspot at cell center (b) Traffic hotspot at cell edge 

Fig. 7: Snapshots of the two traffic hotspots following Log¬ 
normal distribution: The first is close to cell center (a) with 
(ft, a) = (—2,0.5) and the second is at cell edge (b) with 
(l i i cr ) = (—0.75,0.1). 

The parameters ft and a of the Log-normal distribution allow 
to deal with different scenarios of traffic hotspots. Besides, ft 
and o can be tuned to precisely approximate the distribution 
of SINR in each cell of a real network. 

Here, we show the SINR distribution for two cases of user 
locations distributions. The first case is a Log-normal hotspot 
close to cell center with parameters (ft, cr) = (—2,0.5), 
whereas the second one is a Log-normal traffic hotspot located 
at cell edge and having parameters (ft, a) = (—0.75,0.1). 
Both hotspots are depicted in Fig. 7. 

In Fig. 8, we show the CCDF of the SINR for the previously 
cited environments (Fig. 8a for outdoor and Fig. 8b for deep 
indoor environment) considering a hotspot close to the cell 
center (Fig. 7a). Simulation Results confirm once again the 
precision brought by the analytical form of SINR CCDF in 
(43). We notice indeed that numerical and analytical curves 
of SINR CCDF are very close. This consolidates also the 
different observations from the study of uniform location 
distribution, notably regarding the negligible influence of 9 on 
the calculation of the SINR distribution and the importance of 
the penetration margin that deteriorates more or less the SINR 
in deep indoor environment. 

In Fig. 9, we present the SINR CCDF for a Log-normal 
traffic hotspot located at cell edge (see Fig.7b). Even if 9 
influences the ISR value in cell edge (see Fig. 2), the curve of 
SINR CCDF shows to be still insensitive to 9 for traffic hotspot 
at cell border because the simulation result considering 9 is 
almost similar to the analytical expression (43) which does 
not account for 9 variation. All locations in a hotspot at cell 
edge suffer from degraded SINR. We note also that the SINR 
CCDF curve of a hotspot close to cell center is more scattered 
than that of a hotspot at cell edge. This is related to the value 
of <j which is higher for the first hotspot. 

The same previous conclusions regarding the impact of the 
penetration margin in indoor environment are still valid for 
hotspot at cell edge. We could say that the SINR CCDF for 
deep indoor is obtained from that of outdoor environment by 
a translation of almost 3 dB to low values. 

V. Conclusions 

This paper has analytically and fundamentally investigated 
downlink interferences in hexagonal wireless networks. Sev¬ 
eral results have been established for the ISR function / 
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Fig. 8: CCDF of the SINR in a cell having traffic hotspot close 
to its center:(/T = —2 and a = 0.5) for different values of b. 



Fig. 9: CCDF of the SINR in a cell having traffic hotspot at 
cell edge: ( fi = —0.75 and a = 0.1) for different values of b. 


and the SINR distribution. Given a location m = re 10 , we 
have proved that the ISR can be written as an absolutely 
convergent Fourier series on 9 with explicit coefficients. The 
first coefficient, which is the average of the ISR over 9, has 
been explicitly evaluated without any approximation. Besides, 
we have proved for omni-directional network that the ISR is 
a slowly varying function on 9. Its average, as a function of 
r, is enough to capture the essential information of / when 
investigating the SINR distribution. In addition to the provided 
closed formulas of ISR, this paper has also given very simple 
approximations, rigorously proved and shown to be better than 
other approximations in literature. Subsequently, the derived 
expressions have been extended for tri-sectorized hexagonal 
networks and for the use of high frequency reuse patterns. 
Moreover, SINR distribution has been explicitly given for 
every location distribution and entails, of course, the inversion 
of the function g : x r]H Q (x, b) + yoX 2b . Numerical results 

have been carried out for two types of user locations: uniform 
and Log-normal. The exact formula of SINR distribution 
approximates quite well the simulated SINR even when user 
locations form a traffic hotspot at cell edge. 

For future works, we will extend the results of this paper to 
mathematically analyze interferences in perturbed hexagonal 
network model: site location is perturbed by a random variable. 
Apart from the fact that any new result for such a network 
model is scientifically interesting, there are additional motiva¬ 
tions for looking on it: 

• Perturbed hexagonal network model approximates more 
than any other model the geometry of the real network 
[10]. Consequently, performances of perturbed hexagonal 


model would be closer to the performances of a given real 
network, 

• In addition to the parameters of the location random 
variable m, considered in this paper, the perturbed model 
introduces a new parameter, named the average pertur¬ 
bation. SINR or throughput distribution for a given real 
cell can be precisely determined with tuning the average 
perturbation and the first moments of the user location 
variable. 


Appendix A 
Proof of proposition 1 


We first replace in (6) by its expression in (5) 

and use (2), then we can write down fi,k,j(m) = t(to)t(to); 
where 


t(to) = 



0k,j l 3 

Dk,j 


-b 


(47) 


and t(to) is its complex conjugate. 

Since < 1, r(m) is expanded to an absolutely convergent 
series 



T(b + n) 
F(fe)F(l + n) 


(48) 


Arranging the product of the two series r(m)r(m) in order 
to get a Fourier series expansion on 9 and using the definition 
of the Gauss hypergeometric function [18, p. 561], we arrive 
at (8). 


Appendix B 
Proof of lemma 2 

Let z be any complex number such that SR(z) > 1, define 


+oo k— 1 


^ = EE 


(k 2 +j ' 2 -jkf 


(49) 


Observe that oj(z)T(z) is the Mellin transform of the func¬ 
tion -d defined for all complex number t, such that 5R(t) > 0, 
by 


+oo k— 1 

i?(f) = ^2J2e~ t{k2+j2 ~ jk) 

k—1 j =0 


+oo +oo 

EE e -t(k 2 +j 2 +jk) ( 50 ) 

k—1j—0 


Now using the Borwein cubic theta function defined by the 
bilateral sum [32], we have 


E«- 

k,jE Z 


+oo +oo 


= 1 + 2 (t) + 2 EE e 

k—1 j =0 


(51) 


Dividing the last term in the right hand side of (51) into 
two parts and manipulating the indexes of the sum give the 
following identity 

+oo +oo 

EE e - t {k 2 +j 2 -jk) = (52) 

k—1j—0 
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It follows, from the last identity and (51), that 

e~ t{k2+j2+jk) = 1 + M{t) (53) 

k,j GTj 

Using the Lambert Series form of Borwein cubic theta 
function [32], we obtain 

= E ( e t(3fc+l) _ 1 ” e t(3fc+2) _ i ) (54 ^ 

k =0 ' ' 

Finally, applying the Mellin transform to both right and left 
hand sides of the last equation yields 

w(z)r(z) = 3-*c(«) (c{z, 1) -{(*, |)) r(z) (55) 

which completes the proof of lemma 2 . 


Appendix C 

Proof of proposition 2 


Let G s be defined as in (4). Given that G s is a periodic 
function on 0 with period 4^, we assume that it has a 
convergent Fourier series expansion 


+oo 

G 3 (6) = ao + 2 a p cos(3pd) (56) 

p= i 

where 

a p = — f G s (0)cos(3p6)d9, 'ip £ Z 
^ Jo 

Let again m be an arbitrary location in S and recall that 
= ( TO _ S'; k j)/\m — Spk,j |. Using the expression of 
pathloss in (5), we have 


L(m)G s (4n, k ,. 

L i,k,j{m) 


= r 2b E ( 


(m- Sp k ,j) 3p 

{ P | c i26+3p 


(57) 


Thanks to the convergence of the Fourier series (56), the 
coefficient a p can be neglected for p > 2. Consequently, we 
have 

L(m)G a (<f>i t k,j) L(m) 2b 5R [(m - S;, fc ,i) 3 ] 

Li, hJ (m) ~ a °Li, k j(m) + aiT \m - S;, M | 2fc + 3 

(58) 


We now sum all individual ISRs arriving at location m. For 
the second term having the coefficient a lt we limit the sum 
only to the first ring k = 1 because it behaves like an ISR 
with path loss exponent equals 2b + 3. To prove the latter, we 
develop r 2b 5R [(m — Spkj) 3 ] \m — Spk,j\~ 2b ~ 3 as we did in 
proposition 1. It follows that 


G{6 - b) 


E 


L{m) 

Li,k,j{m) 


Gs(4>l : k,j ) 


« aof(m, b) + 2«ir 

Replacing Spi# by 5e dl t 
proof. 


0<i<5,/c>l, 

0<j<fc-l 

26^ 5R \(m - -Si.i.o) 3 ] 
1=0 


(59) 


m-S lX0 \ 2b+3 

and m by re l6 , we complete the 


Appendix D 

Proof of proposition 3 


Let g be the function defined on [0,1) by (41). Using the 
series expansion of Hq(x) given in ( 10 ), it is clear that g 
admits an analytic expansion on x = r/S by 


g(x) = (6 rjui{b) + y 0 )x 2b 1 + E InX 2 


(60) 


where 


7n = 


6i/I'(6 + n) 2 u}(b + n) 


(6 puj{b) + y 0 )T(b) 2 T{n + l) 2 
Let y = g{x), equation (60) can be transformed to 


(61) 


6r)uj(b) + y 0 


X +E 


7 nX 


= x 2 + -j-x 4 + 0(x 6 ) 


(62) 


where O is the Knuth’s notation. The second approximation 
is justified for x sufficiently small, e.g., valid for x < 
Omitting the error terms in equation (62) gives a second-order 
equation, in x 2 , which admits two solutions: 


x 


2 _ 
± “ 


C(y,b ) 2 

l±^/l + lC(y,b) 2 


(63) 


where C(y, b) = ( 6<M Z )+ „„ ) " ■ 

Since x = r/S is a positive real number, only the positive 
solution x 2 + is valid. Now, replacing 71 by its expression 
in (61), we obtain the approximation of x = g~ 1 {y) in 
proposition 3. 
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